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Abstract 

We evaluate the next-to-leading-order QCD corrections of 0{a 2 s GFMf) to the 
Standard- Model £ + £~H, ZZH, and W + W~H couplings in the heavy-top-quark 
limit. Exploiting these results together with knowledge of Ap to the same order, 
we analyze a variety of production and decay processes of low-mass Higgs bosons at 
e + e~ colliders. Specifically, we consider H — ► £ + £~ , H — > £ + £"£' + £'~ , e + e~ — ► ZH, 
Z — > ffH, and e + e~ — ► ffH, with / = £, v. We work in the electroweak on-shell 
scheme formulated with G f and employ both the on-shell and MS definitions of the 
top-quark mass in QCD. As expected, the scheme and scale dependences are greatly 
reduced when the next-to-leading-order corrections are taken into account. In the 
on-shell scheme of top-quark mass renormalization, the C(a^Gi?M 4 2 ) corrections act 
in the same direction as the 0(a s G^M t 2 ) ones and further increase the screening of 
the 0(GfM^) terms. The coefficients of {a s /ir) 2 range from —6.847 for the ZZH 
coupling to —16.201 for the £ + £~H coupling. This is in line with the value —14.594 
recently found for Ap. 
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1 Introduction 



After the celebrated discovery of the top quark the Higgs boson is the last missing 
link in the Standard Model (SM). The detection of this particle and the study of its 
characteristics are among the prime objectives of present and future high-energy colliding- 
beam experiments. Following Bjorken's proposal f2|, the Higgs boson is currently being 
searched for with the CERN Large Electron- Positron Collider (LEP1) and the SLAC 
Linear Collider (SLC) via e + e~ — ► Z — ► ffH. At the present time, the failure of this 
search allows one to rule out the mass range Mh < 64.3 GeV at the 95% confidence 
level ||. The quest for the Higgs boson will be continued with LEP2 by exploiting the 
Higgs-strahlung mechanism e + e _ — ► ZH — > ffH. In next-generation e + e~ linear 

supercolliders (NLC), also e + e~ — > u e v e H via W + W~ fusion and, to a lesser extent, 
e + e~ — > e + e~H via ZZ fusion will provide copious sources of Higgs bosons. 

Once a novel scalar particle is discovered, it will be crucial to decide if it is the very 
Higgs boson of the SM or if it lives in some more extended Higgs sector. To that end, 
precise knowledge of the SM predictions will be mandatory, i.e., quantum corrections 
must be taken into account. The status of the radiative corrections to the production 
and decay processes of the SM Higgs boson has recently been summarized 0. Since 
the top quark is by far the heaviest established elementary particle, with a pole mass 
of M t = (180 ± 12) GeV 0, the leading high-M t terms, of 0(G F M?) } are particularly 
important, and it is desirable to acquire information on their quantumchromodynamical 
(QCD) corrections. During the last year, a number of papers have appeared in which the 
two-loop 0{a s GpM^) corrections to various Higgs-boson production and decay processes 
are presented. These processes include H — > ff, with / ^ b @ and / = b || ]H, Z — > ffH 
and e+e' e+e~ -> v e v e H via W + W~ fusion , gg -»• H 0, , and more 



l~l|l . In this paper, we shall take the next step and tackle with three-loop C^a^G^M 2 ) 
corrections. To keep matters as simple as possible, we shall restrict our considerations 
to light Higgs bosons, with Mh <C M t , and to reactions with colourless particles in the 
initial and final states. Such reactions typically involve the £ + £~H, W + W~H, and ZZH 
couplings together with gauge couplings of the W and Z bosons to leptons. We are thus 
led to incorporate the next-to-leading QCD corrections in the low-M# effective £ + £~H, 
W + W~H, and ZZH interaction Lagrangians. This will be achieved in Section 2. 

Recently, the 0{a 2 s GpM^) correction to Ap has been calculated and found to be size- 
able ||13| . This is relevant for present and future precision tests of the standard electroweak 
theory. It is of great theoretical interest to find out whether the occurrence of significant 
0{a 2 s GpM^) corrections is specific to Ap or whether this is a common feature in the class 
of electroweak observables with a quadratic M t dependence at one loop. In the latter 
case, there must be some underlying principle which is able to explain this phenomenon. 
Our analysis will put us into a position where we can investigate this matter for four 
independent quantities. We shall return to this issue in Section 5. 

The complete evaluation of the one-loop electroweak correction to a process which 
involves more than four external particles is enormously intricate. To our knowledge, the 
literature does not contain a single example of such a calculation. However, the so-called 
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improved Born approximation (IBA) |]T4| allows us to conveniently extract at least the 
dominant fermionic loop corrections. As a by-product of our analysis, we shall illustrate 
the usefulness of the IBA for Higgs-boson production and decay in high-energy e + e~ 
collisions. The appropriate formalism will be developed in Section 3. 

This paper is organized as follows. In Section 2, we shall extend the low- Mh effective 
£ + £~H, W + W~H, and ZZH interaction Lagrangians to 0(a 2 s G F Mf). In the G F for- 
mulation of the electroweak on-shell scheme, knowledge of the QCD-corrected W + W~H 
coupling is sufficient to control the related four- and five-point Higgs-boson production 
and decay processes which emerge by connecting one or both of the W bosons with lepton 
lines, respectively. Contrariwise, the corresponding processes involving a ZZH coupling 
receive additional QCD corrections from the gauge sector, which we shall evaluate by 
invoking the IBA in Section 3. In Section 4, we shall quantitatively analyze the phe- 
nomenological consequences of our results. Section 5 contains our conclusions. 



2 Effective Lagrangians 

Throughout this paper, we shall work in the electroweak on-shell renormalization scheme 
15] , with Gf as a basic parameter, and define (? w = 1 — = M^/Mf |fl6| . In particular, 



this implies that the lowest-order formulae are expressed in terms of Gf, c w , s w , and the 
physical particle masses. The self-energies of the W, Z, and Higgs bosons to 0(oIGfM 2 ) 
for zero external four-momentum squared will be the basic ingredients of our analysis. 
While the results for the W and Z bosons are now well established [O, the Higgs-boson 



self-energy requires a separate analysis, which will be performed here. Our calculation 
will proceed along the lines of Ref. JT3| . We shall employ dimensional regularization in 
n = 4 — 2e space-time dimensions and introduce a 't Hooft mass, /i, to keep the coupling 
constants dimensionless. We shall suppress terms containing 7^ — ln(47r), where 7^ is 
Euler's constant. These terms may be retrieved by substituting /i 2 — > 47re _7E /i 2 . In the 
modified minimal-subtraction (MS) scheme [17|], these terms are subtracted along with 



the poles in e. This is also true for the relation between the MS and pole masses of 
the quarks, so that these terms are also absent when the quark masses are renormalized 
according to the on-shell scheme. Since we wish to extract the leading high-M t terms, 
we may neglect the masses of all virtual particles, except for the top quark. As usual, we 
shall take 75 to be anticommuting for n arbitrary. We shall choose a covariant gauge with 
an arbitrary gauge parameter for the gluon propagator. This will allow us to explicitly 
check that our final results are gauge independent. The requirement that the expressions 
for physical observables be renormalizat ion-group (RG) invariant will serve as a further 
check for our calculation. 

Large intermediate expressions will be treated with the help of FORM 2.0 The 
tadpole integrals which enter the one- and two-loop calculations may be solved straight- 
forwardly, even for arbitrary powers of propagators. The three-loop case is more involved. 
After evaluating the traces, the scalar integrals may be reduced by decomposing the scalar 
products in the numerator into appropriate combinations of the factors in the denomi- 



3 



nator. Subsequently, recurrence relations derived using the integration-by-parts method 
T9| may be applied to reduce any scalar Feynman integral to a small number of so-called 



master diagrams, which remain to be calculated by hand. More technical details may be 
found in Ref. 
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Prior to listing our results, we shall introduce our notation. We take the colour gauge 
group to be SU(iV c ); C F = (N 2 — 1)/(2N C ) and Ca = N c are the Casimir operators of 
its fundamental and adjoint representations, respectively. As is usually done for SU(iV c ), 
we fix the trace normalization of the fundamental representation to be T F = 1/2. In 
our numerical analysis, we set N c = 3. We explicitly include five massless quark flavours 
plus the massive top quark in our calculation, so that we have n/ = 6 active quark 
flavours altogether, i.e., we must not consider «j as a free parameter. We denote the 
QCD renormalization scale by \i. We evaluate the strong coupling constant, a s (fj), at 
next-to-leading order (two loops) in the MS scheme, from 

ft lnln(/x 2 /Af^) 



1 



where is the asymptotic scale parameter appropriate for rif 
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(2) 



are the first two coefficients of the Callan-Symanzik beta function of QCD. We define 
a = Ah = a s (/i)/vr, xt = [G F m 2 t (/j,)/87r 2 V2], X t =_{G F M 2 /8n 2 V2), I = ln[/i 2 /m t 2 (/i)], 
and L = ln(/i 2 /M 2 ), where m t (/i) and M t are the MS and pole masses of the top quark, 
respectively, and G F is Fermi's constant. Using the two-loop relation between m t (M t ) 
and M t [^TJ along with the RG equation for m t (/i), we find 



1 + hC F (-3L - 4) + h 2 C F |l 2 (Jc f - + 
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Riemann's zeta function takes on the values £(2) 
7r 4 /90. The numerical constants [O 

4 „, (-k" 
3, 



vr 2 /6, C(3) w 1.202 057, and C(4) 



So 



9>/3 
-3.027009, 



0.260 434, 



5 4 = 16 Li 4 
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C(4)-4C(2)ln 2 2 + -ln 4 2 



-1.762 800, 



(4) 



4 



where CI2 is Clausen's function and LI4 is the quadrilogarithm, occur in the evaluation 
of the three-loop master diagrams. 

In the following , we shall frequently make use of the QCD expansion of Ap through 
0(a 2 GpM 2 ). For the reader's convenience, we shall list it here for N c = 3 and n/ = 6. 
The MS and on-shell results read HI 



Ap « N c x t 

Ap « N c X t [l - 2.859 912 a - a 2 (5.004 846 L + 14.594 028)] 



15, 



1 + o(2 1 - 0.193 245) + a 2 ( — I 2 + 2.025 330 1 - 3.969 560 



(5) 
(6) 



respectively. 

To start with, we shall construct the Iow-Mh effective 



~H interaction Lagrangian 



through 0(a 2 GpM 2 ). In the following, bare quantities will be labelled with the super- 
script 0. The bare £ + £~H Lagrangian reads 



-m°J°e 



(7) 



where v denotes the Higgs vacuum expectation value. The renormalizations of the lepton 
mass and wave function do not receive corrections in 0(o.™GfM 2 ), where n = 0,1,2, 
so that we may replace m' and £° with their renormahzed counterparts. In the Gf 
formulation of the on-shell scheme, we have 1221 



H 



^ = 2 1 / 4 C# 2 ff(l + * u ), 



with 



n 



WW 



(0) 



M w 



+ 1W0) 



(9) 



Here, Iiwwi.1 2 ) an d ^-hh((1 2 ) are the W- and Higgs-boson self-energies for external mo- 
mentum q, respectively, and the subscript u is to remind us that this term appears as a 
universal building block in the radiative corrections to all production and decay processes 
of the Higgs boson. Consequently, the renormahzed version of Eq. (|7|) reads 



C 



UH 



2 1/i G 1 J 2 m e ££H(l + S u ) 



(10) 



The one-loop expressions for Uww(<l 2 ) and n##(g 2 ) have been presented in Ref. |2c 
The leading-order QCD corrections to Ylwwi.q 2 ) and Uniiiq 2 ) for arbitrary quark masses 
have been found in Refs. 0, 0, respectively. The 0{a s GpM 2 ) term of 5 U has indepen- 
dently been obtained in Ref. |J by using the computational technique outlined above at 
the two-loop level. Here, we shall extend this analysis to 0(a 2 GpM 2 ). The 0(a 2 GVM 2 ) 
term of n^iy(0) may be found in Ref. [ 13f| . The Feynman diagrams pertinent to n##(g 2 ) 



in 0{a s GpM 2 ) come in twenty different topologies. Typical examples are depicted in 
Fig. [I]. We shall renormalize the strong coupling constant and the top-quark mass ac- 
cording to the MS scheme. The appropriate counterterms are listed in Ref. [|2lJ. In this 
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way, we obtain 
n' HH (0)=N c x t !^ + 2l-^ + hCi 
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When we combine Eq. ([□]) with the corresponding expression for Ilvtw(O) P~5j] , the ul- 
traviolet divergences cancel, and we obtain 
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(12) 



With the help of Eq. (|3]), we may eliminate rrit(fi) in favour of M t , which leads to 
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5 U = N c X t I - + hC F (-2C(2) - 3) + h 2 C F 
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6 



(13) 



Equation (13) reproduces the 0{G F M 2 ) and 0(a s G F M 2 ) terms found in Refs. 0, 0, 
respectively. We observe that the new 0(a 2 G F M 2 ) term in Eq. QT3D enhances the QCD 
correction and thus supports the screening of the leading-order M t dependence. The 
choice \i = Mt is singled out, since it eliminates the terms containing L in Eq. (|T5D. 
The nonlogarithmic coefficient of (a s /iT) 2 in Eq. ( |13"D is relatively large; it exceeds the 
corresponding coefficient of Ap in Eq. (^) by approximately 11%. If we consider the ratio 
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of the coefficient of (a s /n) 2 to the one of a s /ir, the difference is even more pronounced; 
the corresponding numbers for Eqs. fll3|) and @ are roughly 9 versus 5. 

A phenomenologically interesting application of Eq. (|T0D is to study the effect of QCD 
corrections on T(H — ► £ + £~). The corrections through 0(o 2 GfM 2 ) to this observable 
may be accommodated by multiplying the Born formula [^] with 

K UH = (1 + 5 U ) 2 

= 1 + 25 UJ (14) 

where we have suppressed terms of 0(G 2 F M^) in the second line. This implies that 5 U is 
gauge independent and RG invariant in these orders. In order to avoid double counting, 



the 0{G F M 2 ) term must once be subtracted when the full one-loop correction [B2[ 



is 



included. A detailed numerical analysis will be presented in Section 4. 

Next, we shall derive the 0(a 2 s G F M 2 ) correction to the \ow-M H effective W + W~H 
interaction Lagrangian. In contrast to the £ + £~H case, we are now faced with the task of 
computing genuine three-point amplitudes at three loops, which, at first sight, appears to 
be enormously hard. Fortunately, in the limit that we are interested in, this problem may 
be reduced to one involving just three-loop two-point diagrams by means of a low-energy 
theorem |25| . Generally speaking, this theorem relates the amplitudes of two processes 
which differ by the insertion of an external Higgs-boson line carrying zero four-momentum. 
It allows us to compute a loop amplitude, M(A — > B + H) , with an external Higgs boson 
which is light compared to the virtual particles by differentiating the respective amplitude 
without that Higgs boson, M.(A — > B), with respect to the virtual-particle masses. More 
precisely 



25 



lim M(A -»■ B + H) = - T t^M(A -> B). 



(15) 



where i runs over all massive virtual particles which are involved in the transition A — > 
B. Here, it is understood that the differential operator does not act on factors of m,; 
appearing in coupling constants, since this would generate tree-level interactions involving 
the Higgs boson that do not exist in the SM. This theorem has variously been applied at 
leading order [||, |25j and has even made its way into standard text books p6| . Special 
care must be exercised if this theorem is to be applied beyond leading order. Then, 
it must be formulated for the bare quantities of the theory, and the renormalization 
must be performed after the left-hand side of Eq. (|T5| ) has been constructed 0. The 
beyond-leading-order version of this theorem || has recently been employed to find the 
0(a s G F M 2 ) corrections to T (h -> bb) §, Y (z -> ffH), and a(e+e" -> ZH) 0. A 
comprehensive review of higher-order applications of this and related low-energy theorems 
may be found in Ref. [11]. An axiomatic formulation of these soft-Higgs theorems has 



recently been introduced in Ref. [27|. 



Proceeding along the lines of Refs. [|1(| [11], we find the bare W + W H interaction 
Lagrangian including its genuine vertex corrections to be 



'W+W-H 



2(M° w ) 2 (W+)\W-»y 
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d(m° t ) 2 (M^) 2 



(16) 
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where it is understood that n^v^(O) is expressed in terms of the bare top-quark mass, 
m f °, while all other quark masses are put to zero. We renormalize the VT-boson mass and 
wave function according to the electroweak on-shell scheme by substituting 



IV J 



M w + 5M W , 

w±{i+8z w y'\ 



(17) 



with the counterterms 



n 



WW 



(0), 



5Z W = -U' ww {0). (18) 

For dimensional reasons, SZyy does not receive corrections in the orders that we are 
interested in. Using Eq. (RF), we thus obtain 
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and the non-universal part herein may be calculated from 
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(19) 
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(21) 



In Ref. ||13|| , Iliyw(O) is expressed in terms of renormalized parameters. Thus, we have to 
undo the top-quark mass renormalization |H| before we can apply Eq. (pT|). Then, after 
evaluating the right-hand side of Eq. (|2l|), we reintroduce the renormalized top-quark 
mass and so obtain a finite result for 5^^, which we combine with 5 U to get Swwh- If 
we define the top-quark mass according to the MS scheme, then the result is 
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The corresponding formula written in terms of M t reads 
5 
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We recover the 0(GpM 2 ) and 0(a s GpM 2 ) terms of Refs. 0, |TTJ, respectively. Similarly 
to Ap and 5 U , the 0{a 2 GpM 2 ) term of Eq. (B3j) supports the screening of the one-loop 



M t dependence by the leading-order QCD correction. Here, the coefficient of (a s /7t) 2 is 
by 26% smaller than in the case of Ap, but it, too, is about five times bigger than the 
coefficient of a s /n. 

From Eq. ([19]) it follows on that T(H — ► PF+W 7- ) receives the correction factor 

Kwwh = 1 + 25wwh- (24) 

Thus, both and Swwh are gauge independent and RG invariant to the orders 

that we are working in. The tree- level formula for T(H — > ly+H^ - ) and its full one- 



loop correction may be found in Ref. \2B\. In order for the Higgs boson to decay into a 
W + W~ pair, it must satisfy Mh > 2My/. On the other hand, the high-Mj approximation 
is based on Mh <C M t . Since these two conditions conflict with each other |IJ, the 
application of Eq. ( |i~9"|) to T(H — > W + W~) is somewhat academic. However, the first 
condition is relaxed to Mh > My/ or removed altogether if one or both of the W bosons 
are allowed to leave their mass shells, respectively. In order to avoid gluon exchange 
between the W + W~H vertex and the external fermions, we restrict our considerations to 
leptonic currents. The resulting class of processes includes H — > (W + )*W~ — > l + i>gW~ , 
H -> W+{W-)* -> W+£-v e , H -> {W + )*(W-y -> £ + v e £'-p ei , as well as e+e" -> 
v e i> e (W + )* (W~)* — > v e VeH via W + W~ fusion. The Born formulae for these 1 — > 3, 
1 — > 4, and 2 — > 3 processes may be found in Refs. |29|, [3^, |31j| , respectively. Since Gp 
is defined through the radiative correction to the muon decay, which is a charged-current 
process, the W-boson propagator does not receive radiative corrections in the orders of 
interest here. Therefore, the correction factors of all these processes coincide with the one 
for T(H -> W + W~). 

Finally, we shall treat the ZZH interaction. The procedure is very similar to the 
W + W~H case. Application of the low-energy theorem (|T5D to the bare Z-boson vacuum 
polarization induced by the top quark yields 

r (mwn TT„„rtvn 

(25) 



C ZZH = {MlfZlZ^'' 



(m° t ) 2 d U zz (0) 
d(m° t ) 2 (Mf)2 
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Again, {M%f = Mf + 5M 2 Z , with 5M Z = n zz (0), and Z° = Together with Eq. (§), 
we then have 

Czzh = 2 1 / 4 G 1 J 2 M 2 Z Z^H(1 + 8 ZZH ), (26) 
where S ZZ h = S u + S^ H , with the non-universal part, 



5, 



ZZH 



1 - 
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0\2 
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d(m^) 2 



n Z z(o) 



(27) 



M) 2 ' 

being separately finite, gauge independent, and RG invariant. 

Starting from the expression for U zz (0) listed in Ref. |13] and repeating the steps of 
the W + W~H analysis, we obtain 



5 zzh = N c x t \-- + hC F 



+ C F -15/ 2 + / 



51- 
12C(2) 



2C(2) + f ) + h 2 C F 



259, , . 593\ 
293 

+ 15C(4) - ^C(3) 



155\ 

IT J 

55 



45 4 + 2D 3 
22 



+ n f 



6 



6 

-Z 2 + / 



73. /n x 613\ 
- T C(2) + 24") 



2435 2 

1053 

~ 2 

65 \ 



449 14 79 



6 



3 

383 



C(2) 



18/ 



35 



S 2 + 2C(4) + — C(3) 
2B A -D 3 + ^S 2 



N c x t 



!)>>+>-!,- 

15, 



1 + a(2 / - 2.017 386) + a 2 — / 2 - 4.815 198 / - 1.086 685 



:28i 



in the MS scheme and 
5 



S 



ZZH 



N r X t 



6 



+ hC F (-2C(2) + 15) + h 2 C F 



449 

2435 2 - — C(3) 
6 



46 



+ C F [AB A + 2D3 - 

+ C A (L (-f C(2) 

/ 113 
+ C(2)(401n2-^ 



1053 



443 



2 * 2 + 2 ^) + 3 



C(3) + C(2) -801n2 



C(2) 
335 

ir 



55 



25 4 - £3 + ^S 2 + 15C(4) - ^C(3) 
4 6 



3697\ / /4 A/oN 



36 



10 



49 

y 



-C(3) - 2C(2) - 9 



1019\ 
~9~) 



115 



-~N c X t [l - 4.684 053 a - a 2 (8.197092 L + 6.846 779)1 
6 



(29) 



in the on-shell scheme. The 0(G F M 2 ) and 0(a s G F M 2 ) terms of Eq. ( P§|) agree with those 
found in Refs. |23|, |10| , respectively. Again, the 0(a 2 G F M 2 ) term of Eq. (|29|) reinforces the 
potential of the QCD corrections to reduce the leading-order M t dependence. Comparing 
<5 zzh with Ap, 5 U , and Swwh, we observe that it has the largest a s /ii coefficient but the 
smallest (a s /ir) 2 coefficient, the ratio of the latter to the former only being about 1.5. 
The (a s /ii) 2 coefficient of S ZZ h is by 53% smaller than the one of Ap in Eq. (0). 
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From Eq. (|26|), we infer that T(H — > ZZ) receives the correction factor 

K ZZH = 1 + 25 ZZH . (30) 

The Born formula for T(H — > ZZ) and its full one-loop correction may be found in 
Ref. |23|. Since the condition 2M Z < M H < M t is likely to be unrealistic |J, the high-M f 



n 

approximation underlying Eq. (^) is of limited usefulness for H — > ZZ. We can evade 
this problem by allowing for one or both of the Z bosons to go off-shell. In addition to 
the information contained in Eq. (|26|) , we then need to account for the corresponding 
corrections arising from the gauge sector. However, in order not to invoke unknown QCD 
corrections, we have to restrict ourselves to the inclusion of lepton lines. The form of 
the additional corrections depends on the considered reaction. It is useful to divide the 
phenomenologically relevant processes into three classes: 

(1) H Z*Z -> ffZ, Z -> Z*H -> ffH, and e+e" -> ZH; 

(2) H -> Z*Z* -> fff'f and e+e" -> -> -> //i? (via Higgs-strahlung) ; 

(3) e+e~ -> e + e-Z*Z* -> e+ e -if (via ZZ fusion). 

Here, / and /' stand for neutrinos and charged leptons. The results for H —* ffZ at tree 
level |H| and at one loop ||, for Z — > //if at tree level |32| and at one loop [Q, for 
e + e~ — > ZiJ at tree level and at one loop |3T| , for H — > fff'f at tree level |j50]| , and for 
e + e^ — > //if at tree level 0] are in the literature. In the next section, we shall discuss 



the corrections to these processes in O(a^GpM^), with n = 0, 1, 2. 



3 Corrections from the gauge sector 

The IBA 10 provides a systematic and convenient method to incorporate the dominant 
corrections of fermionic origin to processes within the gauge sector of the SM. These are 
contained in Ap and Aa = 1 — a/a, which parameterizes the running of the fine-structure 
constant from its value, a, defined in Thomson scattering to its value, a, measured at 
the Z-boson scale. The recipe is as follows. Starting from the Born formula expressed in 
terms of a, c w , s w , and the physical particle masses, one substitutes 

a — > a 



1 - Aa' 

c 2 w ^c 2 w = l-s 2 w = cUl-Ap). (31) 



To eliminate a in favour of Gf, one exploits the relation 

V2 _ 

which correctly accounts for the leading fermionic corrections 
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We shall now employ the IBA to find the additional corrections through 0(oi 2 GfM 2 ) 
to the four- and five-point processes with a ZZH coupling, which we have classified in 
Section 2. We shall always assume that the Born formulae are written in terms of Gp, 
c w , s w , and the physical particle masses. The generic correction factor for class (1) reads 



ra 33 



if) 



ZZH\2 



1- Ap 



1 + 25 



ZZH 



v 2 f + a 2 f 



Qf v f ' 
'v 2 f + a 2 f 



Ap, 



(33) 



where Vf = 2If — As^Qf, Vf = 2If — As^Qf, aj = 2If, Qf is the electric charge of / in 
units of the positron charge, // is the third component of weak isospin of the left-handed 
component of /, and we have omitted terms of 0(G 2 F Mf) in the second line. Similarly, 
the correction factor for class (2) is given by |P 



K. 



(//') 



;i + s ZZH f v 2 + a 2 v 2 , + 4 



vj + a 2 v 2 , + a*, 



1 + 25 ZZH + 2 



1 - 4cf 



Q f v f QyVf, 



v 2 f + a 2 f 



+ 



v 2 fl + a 2 f , 



Ap. 



(34) 



Here and in the following, we neglect interference terms of five-point amplitudes with a 
single fermion trace, since, in the kinematic regime of interest here, these are strongly 
suppressed, by Ty/My, with V = W,Z. Such terms have recently been included in a 
tree-level calculation of T(H -> 2V -> 4/) for M H < M w |5| . 

The correction factor for case (3) is slightly more complicated because the electron and 
positron lines run from the initial state to the final state. Allowing for generic fermion 
flavours, / and /', the Born cross section may be evaluated from 



a(ff - ff'H) 



64^^ 



if + df){v 2 fi + a 2 ,) A ± AvfdfVfiCLfiB 



where 



.4 



dx \ dy 



a(x,y) 



[l + s{y-x)/M 2 f 



(35) 



(36) 



and similarly for B, y'i is the centre-of-mass energy, and the plus/minus sign refers to an 
odd/even number of antifermions in the initial state. The process under case (3), with an 
e + e _ initial state, requires the plus sign. The integrands read 



a{x,y) 
b(x,y) 




ln(l + z) 



x 



+ y 2 \y , 



l + z' 



(37) 
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where z = (y/M 2 z )(s - Mjj/x). The inner integration in Eq. (|36| ) has been carried out 
analytically in Appendix A of Ref. [31]. By means of the IB A, we obtain the correction 
factor for Eq. (|35|) as 



K 



{ff) _ (1 + 5 ZZH f (vj + a)){v), + a),)A ± 4v f a f v r a r B 
Ap) 2 (v 2 + a 2 f )(vj, + a 2 f ,)A ± 4v f a f v r a fl B 
4:C 2 „ ( QfVf QfVf \ 




+ 



l + l/r 



(3? 



where 



±4vfafVpa,fiB 
(v 2 f + a 2 )(v 2 f , + a? fl )A' 



(39) 



We wish to point out that, in the limit r — > 0, coincides with K 2 . Detailed analysis 
reveals that r is quite small in magnitude whenever e + e~ — *> e + e~H via ZZ fusion is 
phenomenologically relevant. In fact, if we consider energies yfs > 150 GeV and demand 
that the total cross section of this process be in excess of 10~ 2 fb" 1 , then we find |r| < 1%. 
This concludes our discussion of the additional QCD corrections to the processes under 
items (l)-(3) originating in the gauge sector. 



4 Numerical results 

We are now in a position to explore the phenomenological implications of our results. We 
shall take the values of our input parameters to be My/ = 80.26 GeV, Mz = 91.1887 GeV 
3§, M t = 180 GeV 0, and af\M z ) = 0.118 0.Q The latter corresponds to af\M t ) = 



0.1071, which entails that = 91 MeV in Eq. (0). If we use the one-loop formula for 

a^\n), i.e., Eq. (P with the second term within the square brackets discarded, 
comes down to 41 MeV. 

Any perturbative calculation to finite order depends on the choice of renormalization 
scheme and, in general, also on one or more renormalization scales. It is generally believed 
that the scheme and scale dependences of a calculation up to a given order indicate the 
size of the unknown higher-order contributions, i.e., they provide us with an estimate of 
the theoretical uncertainty. Of course, the central values and variations of the scales must 
be judiciously chosen in order for this estimate to be meaningful. If the perturbation 
series converges, then the scheme and scale dependences are expected to decrease as the 
respective next order is taken into account. This principle has recently been confirmed 



for Ap [|T3|, 38]. Here, we have the opportunity to carry out similar studies for the three 



additional observables S u , Swwh, and 5zzh- Similarly to Ref. |R|, we have presented our 



results in the on-shell and MS schemes as functions of a single renormalization scale, fx. 



1 Note that this value does not include results from lattice computations. 
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In the MS scheme, one could, in principle, introduce individual renormalization scales for 
the coupling and the mass. For simplicity, we have chosen not to do so. It is natural to 
define the central value of p in such a way that, at this point, the radiative correction 
is devoid of logarithmic terms. This leads us to set p = £M t in the on-shell scheme and 
p = £p t , where p t = m t (p t ), in the MS scheme. We may obtain p t as a closed function of 
M t by iterating Eq. (|3|), with the result that 



1 - AHC F + H 2 C F j-12C(2) + 6 + C F 
+ C A [6C(3) + 8C(2)(-31n2 + l) 



1111 



24 



-12C(3) + 6C(2)(81n2 

71" 
12. 



5) + 



199 



+ n f 



4C(2) + 



where A 



3' 
AH 



-A - 6.458 784 A 



(40) 



a,(M t )/ir. For M t = 180 GeV, Eq. ([p yields p t = 170.5 GeV, in good 



agreement with the exact fix point of Eq. (|3|), which is pt = 170.6 GeV. 



Table 1: Relative deviations (in %) of Ap, 5 U , Swwh, and Szzh from the respective one- 
loop results due to their corrections up to 0(a s ) and 0(a 2 ). The renormalization scale 
dependence is investigated by choosing p = £M t , with £ variable. 





Ap/ApW - 1 [%} 




- 1 [%] 


&WWh/ &WWH ~ 1 [%] 


5zzh/5<£I h - 1 [%} 




0{a s ) 




0(a s ) 


0{o? s ) 


0(a s ) 




0{a s ) 




1/4 


-11.68 


-11.88 


-7.34 


-8.65 


-9.33 


-9.35 


-19.13 


-16.58 


1/2 


-10.63 


-11.72 


-6.68 


-8.34 


-8.49 


-9.24 


-17.40 


-16.83 


1 


-9.75 


-11.44 


-6.12 


-8.01 


-7.78 


-9.04 


-15.96 


-16.76 


2 


-9.00 


-11.11 


-5.66 


-7.67 


-7.19 


-8.79 


-14.74 


-16.51 


4 


-8.36 


-10.74 


-5.26 


-7.35 


-6.68 


-8.51 


-13.70 


-16.15 



Table 2: Relative deviations (in %) of Ap, 5 U , Swwh, and Szzh from the respective one- 
loop results due to their corrections up to 0(a s ) and 0(a 2 s ). The renormalization scale 
dependence is investigated by choosing p = £p t , with £ variable. 





Ap/ApW - 1 [%} 


w - 


- 1 [%} 


Swwh/ <5 


WWH 1 [%] 


w4 x L - 1 [%] 




0{a s ) 




0(a s ) 


0(a 2 ) 


0{a s ) 


0(a 2 ) 


0(a s ) 


0{a 2 s ) 


1/4 


-15.63 


-11.15 


-10.70 


-8.25 


-12.96 


-8.67 


-24.09 


-15.15 


1/2 


-10.89 


-11.55 


-6.87 


-8.22 


-8.71 


-9.09 


-17.77 


-16.56 


1 


-8.96 


-11.19 


-5.63 


-7.79 


-7.15 


-8.86 


-14.69 


-16.51 


2 


-8.64 


-10.88 


-5.82 


-7.56 


-7.11 


-8.71 


-13.49 


-16.16 


4 


-9.24 


-10.85 


-6.81 


-7.70 


-7.92 


-8.84 


-13.41 


-15.93 



In Tables H and H, w e investigate the £ dependence of 5 U , Swwh, and 5zzh and their MS 
counterparts, respectively. For comparison, we also include the results for Ap and Ap. To 
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be specific, we consistently evaluate these quantities to leading and next-to-leading order 
in QCD and study their relative deviations from their respective one-loop values, which 
we denote by the superscript (1), e.g., Ap^ = N c X t , etc. Notice that the on-shell and MS 
results coincide at one loop. In our 0(a s ) analysis, we use the one-loop formula for a s (/x) 
with A^L = 41 MeV and omit the 0(a 2 ) terms in Eqs. @ and (f40|). For the time being, 
let us concentrate on the entries for £ = 1 and assess the effect of the QCD corrections as 
well as their scheme dependence. We observe that, in both schemes, the QCD corrections 
are throughout negative, even for S u and Swwh, where the 0(a s ) and 0(a 2 ) terms are in 
part positive. This is due to the fact that we consistently compute all QCD parameters, 
i.e., a s (/i), m t (fi), and /i t , to the orders under consideration. The reduction in x t , which 
occurs as an overall factor in the MS formulae, happens to overcompensate the positive 
effect of these particular coefficients. Inclusion of the 0(a 2 ) terms in (Ap, 5 U , Swwh, $zzh) 
increases the size of the QCD corrections by (17,31,16,5)%, respectively. In the MS 
case, the increments amount to (25,38,24,12)% of the respective 0(a s ) corrections. As 
might be expected, the scheme dependence of the QCD corrections to this quadruplet of 
quantities is dramatically reduced, by (68,55,71,80)%, as we pass from 0(a s ) to 0(a 2 ). 
Let us now also include the other £ values in our consideration. Within each scheme, we 
determine the scale dependence of the QCD correction to a given quantity by comparing 
its largest and smallest values in the interval 1/4 < £ < 4. As expected, the scale 
dependence is drastically decreased when we take the 0(a 2 ) terms into account, namely 
by (66,37,68,87)% and (89,87,93,86)% in the on-shell and MS schemes, respectively. 
The exceptionally small reduction of the scale dependence in the case of 8 U is due to the 
fact that 5 U has the smallest 0(a s ) term and the largest 0(a 2 ) term of all four on-shell 
quantities. 

Table 3: Coefficients of the correction factors in the form of Eq. ( f41~D for the various 
Higgs-boson decay rates and production cross sections discussed in the text. In the last 
line, x = B/A, where A and B are given by Eq. (|36|), and terms of 0(x 2 ) have been 
neglected. 



K 


Ci 


c 2 


c 3 


Kuh 


7/3 


-1.797 


-16.201 


K WWH 


-5/3 


-2.284 


-10.816 


Kzzh 


-5/3 


-4.684 


-6.847 




-2/3 


-7.420 


4.774 


K? 


-1.272 


-5.249 


-4.445 


Kt ] 


1/3 


6.261 


-53.330 




-0.272 


-14.025 


32.824 


Ki U) 


-0.878 


-6.323 


0.113 




-0.878 - 2.353x 


-6.323 + 9.281 x 


0.113 -39.416x 



In the remainder of this section, we shall stick to the on-shell scheme. In Eqs. (|14"D , 
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(p4|), (pO|), (0), ([54]) , and fl38|), we have presented correction factors for various Higgs- 
boson production cross sections and decay rates in terms of Ap, 5 U , Swwh, and 5zzh- It 
is instructive to cast these correction factors into the generic form 



K = 1 + dApW 



I r | 1 + |aL ) - Cyjr 



(41) 



where Cj (2 = 1,2,3) are numerical coefficients. Notice that we have kept the full p 
dependence in Eq. ([41]). We could have written Eqs. (||), (|i~3|), (|23|), and (p9|) in the same 
way. The fact that the coefficient of aL is universal may be understood by observing that 
K represents a physical observable, which must be RG invariant through the order of our 
calculation, and that, to leading order of QCD, K only implicitly depends on /z, via a. In 
fact, the coefficient of aL is nothing but /3 of Eq. (0). The outcome of this decomposition 
is displayed in Table |3|. In the case of K% , we have treated x = B/A, where A and B are 
defined in Eq. (|36"1), as an additional expansion parameter and discarded terms of 0(x 2 ). 
This is justified because, in practice, \x\ <C 1, e.g., for y/s = 300 GeV and Mh = 100 GeV, 
we find x ~ —5.233 • 10~ 2 . While in the case of the three basic corrections, Kmh, Kwwh, 
and Kzzh, C*2 and C3 are both negative, this is not in general so. In fact, in all composite 
corrections, except for K[ , the 0(a 2 s ) terms partially compensate the 0(a s ) ones. In 
iCj , we even find a counterexample to the heuristic rule [0 that, in the Gp formulation 
of the on-shell scheme, the 0{GpM 2 ) terms get screened by their QCD corrections. In 
the latter case, we also encounter a gigantic value of C 3 . Both features may be ascribed 
to the fact, that, in 0[GpM 2 \ the 6zzh and Ap terms of K^"^ largely cancel. The 
extraordinarily large value of C3 in is also accompanied by a suppression of C%. The 
G\ values of and K% are relatively small, too. We are thus in the fortunate position 
that the leading high-Mj corrections to the 2 — > 3 and 1 — ► 4 processes of Higgs-boson 
production and decay with a ZZH coupling, for which full one-loop calculations have not 
yet been performed, are throughout quite small. Thus, there is hope that the subleading 
fermionic corrections to these processes will not drastically impair the situation. However, 
the IBA does not provide us with any information on the bosonic corrections. 



Table 4: Full one-loop weak corrections (in %) to various Higgs-boson decay rates and 
production cross sections and their 0(GfM 2 ) terms. In the last line, we have used 
v/i = 175 GeV. 



Observable 


M H [GeV] 


0(a) weak 




0{G F M 2 ) [%} 


T{H -> T+T-) 


75 


1.792 


2.369 


T(H -> uUZ) 


105 


1.275 


-0.677 


T(H -> £+£-Z) 


105 


-1.220 


-1.292 


T(Z -> uVH) 


65 


0.024 


-0.677 


T(Z -> e+£~H) 


65 


0.296 


-1.292 


a(e + e~ -> ZH) 


75 


-2.293 


-1.292 



In this context, it is interesting to revisit processes for which the full one-loop weak 
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corrections are known and to investigate in how far the 0(GfM 2 ) terms play a dominant 
role there. Here, we are only interested in reactions which already proceed at tree level. 
Specifically, we shall consider Z -> ffH @ for Mu = 65 GeV, H -> t+t~ [f22| and 
e+e- -> ZH for M# = 75 GeV, and # -> ffZ f§ for = 105 GeV, where 
f — u,£. Our analysis of a(e + e~ — > ZiJ) will refer to LEP2 energy, y/s = 175 GeV. In 
all these cases, the quantumelectrodynamical (QED) and weak corrections are separately 
finite and gauge independent at one loop. In Table f|, we compare the full one-loop weak 
corrections to these processes with their 0{GpM 2 ) terms. In the case of H — > t + t~, 
H — > £ + £~Z, and e + e" — > ZH, the 0{GfM 2 ) terms give a reasonably good account of 
the full corrections, while they come out with the wrong sign in the other cases. However, 
the full calculations for M# = 65 GeV give very small results anyway. On the other hand, 
the 0(GpM%) term for H — > vvZ is suppressed due to a partial cancellation between 
&zzh an d Ap in and cannot be expected to dominate the full correction. Whenever 
the full correction is known, it should be included on the right-hand side of Eq. ([4ID with 
the 0{GfMi) term subtracted. In conclusion, the radiative corrections considered in 
Table £| all appear to be well under control. 



5 Conclusions 

In this paper, we have presented the three- loop O^oIGfM 2 ) corrections to the effective 
Lagrangians for the interactions of light Higgs bosons with pairs of charged leptons, W 
bosons, and Z bosons in the SM. While the demand for corrections in this order is certainly 



more urgent in the gauge sector |T3|, where precision test are presently being carried out, 
our analysis is also interesting from a theoretical point of view, since it allows us to 
recognize a universal pattern. In addition to Ap, we have now three more independent 
observables with quadratic M t dependence at our disposal for which the QCD expansion is 
known up to next-to- leading order, namely S u , Swwh, and Szzh- In the on-shell scheme of 
electroweak and QCD renormalization, these four electroweak parameters exhibit striking 
common properties. In fact, the leading- and next-to-leading-order QCD corrections act 
in the same direction and screen the 0(GfM 2 ) terms. Even the sets of a s /7r and (a s /7i) 2 
coefficients each lie in the same ball park. For the choice \x = M t , the coefficients of a s /ir 
range between —1.797 and —4.684, and those of (a s /ir) 2 between —6.847 and —16.201. If 
we compare this with the corresponding coefficients of the ratio fif/M 2 , which are —2.667 
and —11.140, then it becomes apparent that the use of the top-quark pole mass is the 
origin of these similarities. Here, fit = m t (fit), for which we have presented a closed 
two-loop formula. If we express the QCD expansions in terms of fit rather than M t and 
choose ft = fi t , then the coefficients of a s /7r and (a s /ii) 2 nicely group themselves around 
zero; they range from —2.017 to 0.870 and from —3.970 to 1.344, respectively. This 
indicates that the perturbation expansions converge more rapidly if we renormalize the 
top-quark mass according to the MS scheme. Without going into details, we would like 
to mention that the study of renormalons [091 offers a possible theoretical explanation of 



this observation. Since the on-shell and MS results coincide in lowest order, this does, of 
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course, not imply that the QCD corrections are any smaller in the MS scheme. It just 
means that, as a rule, the O(GpM^) terms with M t replaced by the two-loop expression for 
fit are likely to provide fair approximations for the full three-loop results. Furthermore, we 
have demonstrated that, similarly to Ap, the scheme and scale dependences of 5 U , Swwh, 
and Szzh are considerably reduced when the next-to-leading-order QCD corrections are 
taken into account. Armed with this information, we have made rather precise predictions 
for a variety of production and decay processes of low-mass Higgs bosons at present and 
future e + e~ colliders. In all the cases considered here, the radiative corrections appear to 
be well under control now. 
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